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Dissecting the role of initial collision geometry for jet quenching observables in 
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The observation of large azimuthal anisotropy or V2 for hadrons above pT > 5 GeV/c in Au+Au 
collisions at ^ysnn = 200 GeV has been a longstanding challenge for jet quenching models based on 
perturbative QCD (pQCD). Using a simple jet absorption model, we seek to clarify the situation 
by exploring in detail how the calculated V2 varies with choices of the collision geometry as well as 
choices of the path length dependence and thermalization time tq in the energy loss formula. Besides 
the change of eccentricity due to distortion from gluon saturation or event-by-event fluctuation, 
we find that the V2 is also sensitive to the centrality dependence of multiplicity and the relative 
size between the matter profile and the jet profile. We find that the V2 calculated for the naive 
quadratic path length dependence of energy loss, even including eccentricity fiuctuation and the 
gluon saturation, is not enough to describe the experimental data at high pr (~ 6 GeV/c) in 
Au+Au collisions. However, it can match the full centrality dependence of «2 data if higher power 
path length dependence of energy loss is allowed. We also find that the calculated V2 is sensitive 
to the assumption of the early time dynamics but generally increases with tq, opposite to what 
one expects for elliptic fiow. This study attests to the importance of confining the initial geometry, 
possibly by combining jet quenching V2 with elliptic fiow and other jet quenching observables, for 
proper interpretation of the experimental data. 

PACS numbers: 25.75.-q 
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I. INTRODUCTION 

After the discovery of strongly interacting Quark 
Gluon Plasma (sQGP) at the Relativistic Heavy Ion Col- 
lider (RHIC) in 2005 [1], the focus of the heavy ion 
community shifted toward a detailed characterization of 
the properties of the sQGP. One of the primary tools 
is jet quenching or the suppression of high transverse 
momentum (px) hadron yields as a result of in- medium 
radiative energy loss of high px jets 0t3|- Due to the 
large momentum scale of the jets and asymptotic free- 
dom of Quantum Chromodynamics (QCD), jet quench- 
ing is usually thought to be described by the pertubative 
QCD (pQCD) framework, which assumes that jets cou- 
ple weakly with the medium, even though the medium 
itself is strongly coupled. Jet quenching models based on 
pQCD have been developed to describe measurements 
on single hadron yield [3, M, di-hadron correlation [3, Q , 
and 7- hadron correlation |9l4l2l| . Initial estimates of the 
properties of sQGP, such as the momentum broadening 
per mean free path, q = (fc|.)/A, and energy loss per unit 
length, dE/dl, have been obtained |13| . 

Despite its early successes, the pQCD description of 
jet quenching faces several challenges (see Ref [ij]). One 
observable that has thus far defied the pQCD descrip- 
tion is high pt V2 or azimuthal anisotropy of particles 
emitted relative to the reaction plane (RP) in Au-|-Au 
collisions, dN/d{4> - *rp) oc (1 + 2w2Cos2((/) - *rp)). 
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FIG. 1: Figure adapted ft-om Ref. [jj]. Data points: PHENIX 
tt" V2 at Pt > 6 GeV/c; open symbols: three pQCD model 
calculations taken from Ref. [l3 |: lines: geometric model cal- 
culations with different assumptions on path-length depen- 
dence M,^. 



Such azimuthal anisotropy ensues because the jet yield 
is more suppressed along the long axis of the fireball 
(out-of-plane) than the short axis (in-plane). Thus, the 
V2 value is sensitive to the path length (I) dependence 
of energy loss, which scales, in the pQCD framework, 
as /S.E ex I and IS.E (x P for elastic and Landau- 
Pomeanchuk-Migdal(LPM) radiative energy loss [l^, re- 
spectively. Currently, most pQCD models undershoot 
the V2 value by as much as factor of 2 in the experi- 
mentally accessible pr range (pr < 10 GeV/c) [la,ll3]- 
We illustrate this situation with Fig. [1] borrowed from 



Ref. [l8|, which compares three mainstream pQCD model 
calculations (abbreviated as AMY, ASW, and HT) ^ 
with recent precision PHENIX data at pT ~ 6 GeV/c. 

This, together with its failure in describing heavy fla- 
vor suppression [19|, call into question the perturbative 
assumption used in the pQCD framework. It may hap- 
pen that the coupling between the jet and the medium 
for typical RHIC jet energy of pr < 10-20GeV/c is 
still strong enough [23, such that path length depen- 
dence and the color charge dependence are modified from 
pQCD expectation. In fact, calculation based on anti-de 
Sitter/conformal field theory (AdS/CFT) technique for 
strongly coupled plasma suggests that AE ex P [2l|, [22| 
and q ex y/asYuNc Q: instead of AE ex P and q ex 
asN^ for pQCD. This higher order path length depen- 
dence could explain the large anisotropy ^24]. Liao and 
Shuryak [23 argue that most energy loss in sQGP is con- 
centrated around Tc] such a non-monotonic dependence 
of energy loss with energy density apparently achieves 
better description of the data, as shown by Fig. [TJ 

It is tempting to conclude from this discussion that 
the data favor a I dependence stronger than the naive 
AE ex l"^ implied by the pQCD radiative energy loss. 
However, as was pointed out in Ref. [l3|, the magnitude 
of the anisotropy is also very sensitive to the choice of 
initial collision geometry, which is poorly constrained. 
The collision geometry used by most jet quenching cal- 
culations is obtained from the so-called Optical Glauber 
model [26|, which assumes a smooth Woods-Saxon nu- 
clear geometry for Au ions. It ignores two important 
modifications: an event-by-event distortion of the shape 
of the overlap from random fluctuation of positions of 
participating nucleons [28]; and a possible overall distor- 
tion of the shape of the overlap due to, e.g. gluon satu- 
ration effect (so called CGC geometry [23]). Both effects 
are shown to lead to 15%-30% corrections in the hydro- 
dynamic calculation of elliptic flow at low pT [30| ; they 
were also shown in Refs. [23,|31| to play an important role 
for jet quenching calculation of azimuthal anisotropy at 
highpT- 

Furthermore, the way that collision geometry influence 
the jet quenching V2 is quite different from that for hydro- 
dynamic description of low px V2 ■ Hydrodynamic flow is 
a self generating process driven by the shape or eccentric- 
ity (e = ; 2\ I / 2\ ) of a single matter profile, i.e. V2 ex e; 
whereas the V2 from jet quenching requires both the pro- 
file for the bulk matter AND the jet production points. 
The two profiles may not necessarily have the same spa- 
tial distribution because various nuclear effects at initial 
state may induce sizable momentum (e.g., Bjorken mo- 
mentum fraction 2pt/\/s) and position dependent modi- 
fication, analogous to the generalized parton distribution 
for proton. Hence high pT V2 depends not only on the 
eccentricity of the fireball, but also on the matching (rel- 
ative size and shape) between the jet and the matter 
profiles. Understanding the role of geometry and scal- 
ing behavior of the data such as those in Ref. [32, |33| is 
important for proper interpretation of the experimental 



data. 

In this article, we investigate the sensitivity of the jet 
quenching V2 on the choices and uncertainties of the col- 
lision geometry for the bulk matter. We check explicitly 
the scaling and violation thereof with the bulk eccen- 
tricity. We explore, in the context of these uncertain- 
ties, whether the data allow for high order I dependence 
of energy loss. The prospects of constraining the initial 
collision geometry using V2 and other jet quenching ob- 
servables, such as single inclusive suppression Raa, inclu- 
sive away-side suppression /aa and associated anisotropy 
Wn**, are discussed. 



II. MODEL IMPLEMENTATION 

We generate the Glauber geometry using an improved 
version of the publicly available PHOBOS code 34]. 
Each Au ion is populated randomly with nucleons with 
a hard-core of 0.3 fm in radii, according to the Woods- 
Saxon distribution with a radius of 6.38 fm and diffuse- 
ness of 0.535 fm. A nucleon-nucleon collision is consid- 
ered to happen when their distance in the xy-pleme fall 
within x/o'n'JfV''' = 1-16 fm (hard-sphere assumption), 
corresponding a n-n cross-section of o',\"{''=42 mb. Sub- 
sequently, the number density of nucleons participating 
in the collision (ppart(a;,y, &)) and the number density of 
binary collisions {pcoiiix,y,b)) can be determined in the 
xy-plane as function of impact parameter b. Here the x 
direction is always chosen to be along the line connecting 
the centers of the two ions. Denoting Ta as the thickness 
function for Au ion, they can be approximated with the 
following expression when nucleon size is ignored. 



Ppiirt{x,y,b) 



Pcoiiix, y,b) 



TA[x + -,y 



l-P{x+-,y 



l-P{x+-,y 



Ta[x- -,y 

ai:^,^'TA[x + ^,y^TA(x-^,y],{l) 



where P{x, y) = [1 



'Ta{x,v) 



and A = 197 is the 



number of nucleons in Au ion. 

We generate the CGC geom etry using the MC-KLN 
model by Drescher & Nara |27l . |29| , which is based on the 
well known KLN (Kharzeev-Levin-Nardi) kx factoriza- 
tion approach [3^ . In a nutshell, the MC-KLN model cal- 
culates the CGC geometry event by event by modifying 
the output from a Monte-Carlo Glauber model. Specifi- 
cally, the transverse gluon density profile, dn/dy(x,y,b), 
is calculated through the kx factorization formula, with 
the saturation scale Q^ of each Au ion set to be propor- 
tional to its thickness function Ta or Tb- To ensure inter- 
nal consistency, the MC-KLN code is adapted to the same 
Glauber algorithms as the PHOBOS code (same hard- 
core nucleons and identical Woods-Saxon parameters). 
The obtained gluon density scales approximately [231 as 



inin{TA, Tb} in the x direction and 1/2{Ta + Tb) in the y 
direction, which leads to a 20%-30% increase of the eccen- 
tricity relative to the Glauber geometry (see left panels 
of Fig. [S]). When implemented in hydrodynamic model 
calculations [30, |36[ , a similar amount of increase is seen 
for the predicted elliptic flow signal. 

We account for initial geometry fluctuation in Glauber 
geometry by re-centering and rotating all participants, 
such that the "participant plane" (PP), defined as the 
minor axis direction of all participants (see Ref. |28|). 
aligns with the lab frame. We then sum all events to- 
gether to give the overall participant density profile. The 
same amount of shift and rotation is then applied for 
all binary collisions to get the overall density profile for 
jet production points. We repeat the same procedure 
to the gluon density profile for CGC geometry. Impor- 
tant variables include the orientation of the participant 
plane (^part) for either participants or gluon density ^, 
eccentricity with respect to the reaction plane (crp), ec- 
centricity with respect to the participant plane (cpart), 
and average root mean square (RMS) size of the ellip- 
soid ((Tt-)- They are calculated for each event as 



tan(2* 



part^ 
^part 



^^ly 



K-^D' + 



4cr2 
^^ xy 



- <^x = ^y + <^x I 






(2) 



a^ and a^y are the event-by-event 



(co)variances of participant density profile for Glauber 
geometry or gluon density profile for CGC geometry, re- 
spectively, and o'^ and o'y are variances defined in the 
rotated frame. We emphasize that the participant plane 
angle, ^part, should be the natural frame for both hy- 
drodynamic flow and jet quenching. However, it is tilted 
by a different amount in the case of the CGC geometry 
from the Glauber geometry. 

In this work, the magnitude of the jet quenching 112 
depends on the following four control factors: 

1. Energy loss formula, including the path length 
dependence, thermalization time etc. 

2. Eccentricity, including event-by-event fluctuation 



and the shape of collision geometry (e.g. CGC vs 
Glauber) . 

3. Centrality dependence of the total multiplic- 
ity. Because the jet quenching strength is fixed in 
most central collisions, if matter density falls faster 
toward peripheral collisions, we expect less suppres- 
sion and smaller vi in peripheral collisions. 

4. The size of the matter profile relative to the 
jet profile. If the transverse size of the matter 
profile is smaller than that for the jet profile, more 
surviving jets should originate from the corona re- 
gion, leading to a smaller vi. 

Clearly, the collision geometry (items 2-4) plays an es- 
sential role for proper understanding of the energy loss 
mechanism (item 1). In contrast to hydrodynamic de- 
scription of low px "Vi , which depends only on the eccen- 
tricity of the ellipsoid, jet quenching description of high 
•pT V2 is sensitive to two more aspects of the collision ge- 
ometry (items 3 and 4). The primary goal of this work is 
to understand how the jet quenching V2 depends on the 
underlying choices of eccentricity, centrality dependence 
of multiplicity, and matching between the matter and the 
jet profile. 

We base the study on the following four matter pro- 
files (three versions of Glauber geometry and one CGC 
geometry) 

Po{x,y,b) = ppart(a;,y,6), 
pi{x,y,b) = pcoi\{x,y,b), 

1 A 

P2{x,y,b) = -^—Ppart{x,y,b) + 6pcoii(x,y,b), 
P3{x,y,b) = pcGc{x,y,b) ^ dn/dy{x,y,b), 



with the corresponding integral form 



dxdy poix,y,b) 

dxdy pi{x,y,b) 

dxdy p2{x,y,b) 

dxdy p3{x,y,b) 



Np^rtib), 



Ncouib), 
l-S 



Np^rtib) + SN^oiiib) 



dN/dyib), 
dN/dyib), 



(3) 



^ We emphasize that the PP (and 'tpart) is calculated in the posi- 
tion space of nucleons or gluons in a simulated collision, whereas 
experiments measure the so called event plane (EP) using final 
state particles in momentum space of a real collision. Both PP 
and EP include fluctuations and approximate the true RP of 
their respective collisions, but they may not coincide with each 
other. 



where Ppart and pcoU are transverse participant density 
and collision density from the Glauber model, respec- 
tively: P2 is the two component Glauber model from 
Ref. [231 with (5 = 0.14 [30], and ps = pcGC is the trans- 
verse gluon density from MC-KLN. Both p2 and ps have 
been adjusted [20 such that their total integrals match 
the centrality dependence of the charged hadron multi- 
plicity, dN/dy{b), at RHIC ^. For the first two profiles. 



po and pi , we can enforce the same centrality dependence 
as dN/dy{b) by applying a centrality dependent scale fac- 
tor: 

Mui. ,. dN/dy{b) 

Po ix,y,o) = — — T7-ppart(a;,y,b), 

J Vpart (0) 

dxdy pT\x,y,b) - dN/dy{b), 



dxdy pf''\x,y,b) - dN/dy{b), 



(4) 



These scale factors essentially account for the different 
centrality dependence trends between iVpart and iVcoii rel- 
ative to dN/dy. Figure [5] shows the rate of change of 
dN/dy and iVcoii relative to iVpart normalized to unity 
for most central points. Clearly, the A'coii has the fastest 
change vs. centrality, followed by dN/dy, and iVpart has 
the slowest change vs. centrality. Nevertheless, the re- 
sulting profiles, p^^^ and pj^"', still maintain their origi- 
nal shape and size. They are used to study the sensitivity 
of V2 to the centrality dependence of the multiplicity. 
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FIG. 2: (Color online) The centrality dependence of 
dN/dy/Nps,rt (filled) and A^coii/A^part (open) normalized by 
their values in most central 0%-5% bin. This plot illustrates 
different rates of change as a function of centrality, between 
A''part, A^'coii, and dN/dy. 



participant planes ^, 

P^°*(X, y) = Pn{x cos *part - V sin *part, 
X sin *part + y COS *part), 

or by readjusting the overall size by a constant scale fac- 
tor a to match to that of po (see discussions in Secs. lIIIBJ 
and ImC)) 

p^'='''"'=(x,j/)=p„(aa;,ay), 

or by a combination of matching multiplicity, rotation or 
readjusting size. In the case that a profile is obtained 
via several operations, we use appropriate superscript to 
indicate that. For example, Pq° ' " indicates the mat- 
ter profile obtained by rotating the event-by-event par- 
ticipant profile according to the participant plane angle, 
followed by matching its total integral to dN/dy for each 
centrality bin. Note the order of these operations has no 
significance because they factorize. 

We implement jet quenching using the simple jet ab- 
sorption model of Ref. ilTi] . It provides a transparent way 
of investigating the sensitivity of jet quenching observ- 
ables to choices of the collision geometry. In this model, 
back-to-back jet pairs are generated according to the bi- 
nary collision density profile in xy-plane with uniform 
orientation. These jets are then propagated through the 
medium whose density is given by matter profile p(x, y), 
with a survival probability e~'^^ . In the default setup, 
matter integral / is calculated as 



/ = 



dl 



I 



l + li 



-p{-^ + {l + lo)^) 



dl p(^ + ;v)(5) 



for a jet generated at 1^ = {x, y) and propagated along 
direction v. This corresponds to a quadratic depen- 
dence of absorption (ex Idl) in a longitudinal expanding 
or 1-l-lD medium (ex 1/{Iq -f I)) with a thermalization 
time of la ~ ctq. It is fixed to by default but we ex- 
plore non-zero value of ^o in Sec. IIIIEI The absorption 
coefficient k (which controls the jet quenching strength) 
is chosen to reproduce Raa — (e^"^) ^ 0.18 for 0%- 
5% 7r° data [4l|. We explore path length dependence by 
extending Eq. [S]to four different functional forms. 



/„. = / dl r-^ 

h 



p(^ + Zv), m = 1,2,3,4, (6) 



where m = 1 and m = 2 correspond to I dependence for 
radiative, and AdS/CFT energy loss in 1+lD medium, 
respectively. 

A typical calculation starts by choosing one of the four 



Besides the two profiles obtained by matching to 
dN/dy, we are also interested in several other variants 
of Po — P3, obtained either by rotation of their respective 



^ Note that the participant plane angle 'I'part is defined separately 
for po-P3 • For Pi , it is determined by the minor axis of all binary 
collisions, while that for p2 is determined by both participants 
and binary collision with appropriate weights given in Eq. [3] 



matter profiles (Glauber geometry po — p2 or CGC geom- 
etry Pa) and applying appropriate modifications (specify- 
ing Rot, Mul, and/or Resize). We specify the jet absorp- 
tion scheme by varying thermalization time /q or the or- 
der of path length dependence m. We then fix the k, value 
by matching i?AA ^ 0.18 in most central collision (We ex- 
plore the uncertainties of k arising from experimental un- 
certainties of i?AA, and discuss their implications in Ap- 
pendix [B]). On the other hand, the jet production profile 
is always given by Pcoii- We stress that k is the only free 
parameter, and has similar role as the g, and it is tuned 
independently for each one of these running modes (there 
are ~ 100 of them, depending on the choice of matter 
profiles, to, and modifications of matter profiles). Once 
K, is known, we can predict the centrality dependence of 
the single hadron suppression (i?AA), jet quenching V2 
which can be expressed as V2 — (e^"^ cos 2 (0 — ^part)); 
and away-side per-trigger yield suppression (/aa)- 

Finally, we point out that the U2-like modulation is 
found to be the dominating contribution to the azimuthal 
anisotropy obtained in our calculations. The higher or- 
der terms, mostly U4, are found to be less than 10% of U2 
value for all running modes. Thus we can safely assume 
that the azimuthal distribution of particle production rel- 



ative to the PP angle follows a 1 
shape (for example see Fig. [T9)) . 
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III. RESULTS 
A. Glauber geometry based on participant profile 

As mentioned previously, this work investigates three 
versions of Glauber geometry, participant profile pq , col- 
lisional profile pi and two component profile p2 and their 
variants, po is our default Glauber geometry and is the 
topic of this section; we shall discuss pi and p2 in Sec- 
tion HirOl 

Figure [3] shows the V2 calculated for Glauber geom- 
etry p^"\ that is, participant profile scaled to match 
the experimental multiplicity. Results are presented in 
left panel for four different path length dependencies 
{I1—I4 from bottom to top) with (solid lines) and without 
(dashed lines) taking into account the fiuctuation of PP 
angle. They are compared with the PHENIX tt" V2 data 
integrated above 6 GeV/c from Fig[T] The right panel 
shows the ratios of calculated V2 for /i — I4 (solid lines) 
and the ratio of the eccentricity (open circles) between 
with and without including the fiuctuations. 

We see that increasing m (the order of I dependence) 
significantly increases the V2 for mid-central collisions, 
but they all systematically under-predict the data to- 
ward central collisions. In fact, the calculated V2 for cen- 
tral collision is insensitive to the functional form of path 
length dependence, due to the almost isotropic shape 
of the overlap. This situation is dramatically improved 
when the fluctuation in the PP angle is included. The rel- 
ative increase in V2 is about 15% for mid-centrality, and 







1 1 1 1 1 1 1 

pMul 




0.2 


: — 


qRoI.MuI 


'■ 


0.15 


■ •</- 




I4' 


> 




'A 


::>^ 


! ^ 


0.1 


- /-y 






- 






~-\X% 


^ ^' ■ 


0.05 




1 1 1 1 1 1 1 


*^," 




100 200 300 

Npar, 



100 200 300 

Npar. 



FIG. 3: (Color online) Left: V2 calculated with po""' (par- 
ticipant density profile scaled to match dN/dy) with (solid 
lines) and without (dashed lines) rotation to the participant 
plane compared with data at 6 GeV/c (solid circles) for I1—I4, 
(from bottom up); Right: Corresponding ratios between with 
and without rotation for t;2S from /i — I4 (lines) and for the 
eccentricity (open circles). 



is signiflcantly larger for central and peripheral collisions. 
This is consistent with previous studies of low px V2 or 
elliptic flow, which shows that PP fluctuation needs to be 
included in Gu+Cu and central Au-|-Au collisions in or- 
der for hydrodynamic model prediction to work [28l . |42| | . 
It is interesting to see that the fractional increase of V2 
for /i is similar to the fractional increase in eccentricity 
(i.e., epart/eRP in the right panel). However, the ratios 
indicate that the fractional increase for iVpart > 50 is 
successively larger for larger m. This is because larger 
m places more weight to the large I region; and thus is 
more sensitive to changes in shape. 
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FIG. 4: (Color online) Left: «2 calculated with po°' (partici- 
pant density profile rotated to participant plane) with (solid 
lines) and without (dashed lines) scaling to match dN/dy for 
I\ — li (from bottom up); Right: Corresponding ratios be- 
tween with to without scaling to dN/dy for W2S from /i — 1 4 
(lines) and for the eccentricity (open circles). Note that the 
eccentricities are identical for the two cases. 

The high pT azimuthal anisotropy, being the result of 



jet quenching, depends not only on the shape, but also 
on the average density or total multiplicity of the matter 
profile. To illustrate this point. Fig. |4] shows the V2 cal- 
culated for matter density /5q-°* and Pq° ' " . They have 
identical shape and size for each centrality selection, but 
the integral of Pp °*' " drops more rapidly to lower iVpart . 
Because jet absorption strength k is tuned to reproduce a 
common suppression level in central collision, the profile 
whose average density varies more rapidly with centrality 
is expected to show less suppression and less V2 in periph- 
eral collisions. Indeed, the calculated V2 for Pq° ' " is 
smaller than that for Pq °' due to a faster fall off toward 
peripheral bin. 

Figure 2] also shows a weakening of the sensitivity for 
larger m. This is because the weighting from the large I 
region is reduced in peripheral collisions due to a smaller 
geometrical size. That reduction is stronger for larger m, 
which leads to smaller sensitivity for large m. 

One may argue that since dN/dy{b) is constrained by 
experimental data, there should be no uncertainty asso- 
ciated with the modeling of centrality dependence. How- 
ever, the matter profiles that were tuned to Au-|-Au 200 
GeV data typically shows ~ 10% deviation from Cu-|-Cu 
or Au-(-Au at different collision energies [4j, |4J|. Fur- 
thermore, many current pQCD model calculations use 
profiles that do not match the dN/dy data. For example, 
various 1-l-lD energy loss models assumes energy loss or 
q to be proportional to either ppart [8, 45, 461 or PcoU [43]; 
Recent more sophisticated calculations [iljl, [ij, |4^ based 
on 3D-I-1 hydrodynamics model of Nonaka and Bass p9j . 
assume the energy loss or q oc C''^ with e ex O.Gpcoii + 
0.4ppart [43 J which is also different from dN/dy. So it 
seems reasonable to use the difference of the V2 in right 
panel of Fig. |3] as one of the uncertainties in theoretical 
implementation of initial geometry. 

Figure [3] and |4] represent the general style of the pre- 
sentation of the V2 calculation in this article: The left 
panels always shows the V2 values for /i — /4 compared 
between two matter profiles with (solid lines) and with- 
out (dashed lines) a particular geometrical effect; the 
right panels always show the ratios between the two (solid 
line divided by dashed line). In most cases, the ratio of 
their eccentricities is shown as open circles on the right 
panel to compare with V2 ratios. Finally, we stress that 
the jet production profiles are always sampled from pcoii 
throughout this study, so any difference in the calculated 
V2 can be attributed to the differences between the two 
matter profiles. 



B. CGC geometry 

As outlined in the introduction, the MC-KLN model 
is built on the standard Monte-Carlo Glauber model. So 
CGC matter profile ^3 contains both the overall modifi- 
cation of shape due to gluon saturation, and the event- 
by-event fluctuation stemming from participant fluctua- 
tion at Glauber level. In addition, we can safely use the 



same binary collision profile for the jet production, given 
that the saturation effects are not expected to modify 
hard processes with momentum transfer well above the 
saturation scale, Q^ ^ Ql- 

Another important feature of the CGC geometry via 
the MC-KLN model is that despite having a larger ec- 
centricity, its overall size is about 4%-8% smaller than 
the participant profile. One can see it quantitatively in 
Fig. [5l which compare the eccentricity and overall RMS 
width ((Tj.) between p^ and pg. This narrowing of CGC 
geometry was pointed out before by the authors of MC- 
KLN model (see the preprint version of [27[), and can 
be seen more clearly by plotting the ID projections of 
medium profiles along the x (in-plane) and y (out-of- 
plane) directions (Fig. ^. The projections show that 
CGC profile is narrower than Glauber profile in both x 
and y direction, however since a^ is reduced more than 
ay in MC-KLN vs Glauber, epart is larger in MC-KLN. 
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FIG. 5: (Color online) Left: Eccentricities for Glauber geome- 
try calculated from participant density profile po (solid circles) 
and for CGC geometry pa (open circles) in the top panel and 
the corresponding ratio in the bottom panel. Right: Same as 
right panels except they are for the RMS size ar 
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FIG. 6: (Color online) The Glauber geometry (po) and CGC 
geometry (pa) projected onto the x and y axes for 15%-20% 
and 30%-35% centrality bins. They are normalized to 1 at 
the maximum. 

Figure [7] compares the V2 calculated for CGC geome- 



try (pa) and Glauber geometry (po) in their respective 
rotated frames. The CGC geometry does lead to a larger 
V2', however, the amount of increase is only half of the 
increase in eccentricity. To check whether the breaking 
of the eccentricity scaling can be attributed to the 4%- 
8% mismatch between the two profiles, we re-scale the 
RMS size of the CGC geometry to match that for the 
Glauber geometry for each centrality while preserving its 
original shape. Figure [8] shows the same comparison af- 
ter the scaling is applied and k is re-tuned. The ratios 
of the calculated V2 now match well with the ratio of the 
eccentricities. 
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FIG. 7: (Color online) Left: V2 calculated for Glauber ge- 
ometry (po^°') and CGC geometry (pf^"*) in their respective 
rotated frames for Ji — I4 (from bottom up). Right: Corre- 
sponding ratios for «2S from Ji — I4, (lines) and for the eccen- 
tricity (open circles). Note that the CGC geometry has larger 
eccentricity. 
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FIG. 8: (Color online) Same as Fig. [3 except that the RMS 
size of CGC geometry has been stretched to match Glauber 
geometry po (according to ratio shown in bottom right panel 
of Fig. 0) 



These comparison plots clearly show that both eccen- 
tricity and the size contribute to the difference of the V2 
between CGC and Glauber geometry: While the eccen- 
tricity of the CGC geometry increase by about 10%-30% 



relative to Glauber geometry, its transverse size shrinks. 
The latter change increases the fractional contribution of 
surface jets, which have smaller V2- In contrast, there is 
no such bias for hydrodynamic calculation of low px V2, 
which depends on the shape, not the size, of the matter 
profile. 



C. Glauber geometry based on collision profile and 
two component profile 

The preceding discussion alludes to an interesting pos- 
sibility: For a given energy loss formula and jet produc- 
tion profile, as long as the matter profile is adjusted to 
a common reference ar and dN/dy, the V2 depends only 
on the eccentricity of the matter profile. Here we further 
test this ansatz by using a matter profile pi that is very 
different from pg. Comparing to po, pi has much larger 
eccentricity (left panels of Fig. [S]) which should increase 
the calculated V2- On the one hand, it has stronger cen- 
trality dependence of total integral (Fig. [2]) and 10%-15% 
smaller ar (right panels of Fig. ^ , both are expected to 
significantly decrease the calculated V2- 
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FIG. 9: (Color online) Left: Eccentricities for Glauber geom- 
etry calculated from participant density profile po (solid cir- 
cles) and from collision density profile pi (open circles) in the 
top panel and the corresponding ratio in the bottom panel. 
Right: Same as right panels except they are for the RMS size 
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Figure [TUl compares the V2 calculated for pi and po in 
their respective rotated frames. The calculated V2 falls 
well below the experimental V2 data in central collisions, 
which suggests that the Glauber geometry based solely 
on collision density profile with eccentricity fluctuation is 
ruled out. The decrease of the V2 is largely attributable to 
the stronger centrality dependence and smaller size of pi , 
and can be seen more quantitatively in the right panel, 
which appears as a large suppression of the f 2 ratios from 
the expected eccentricity ratio. 

To dissect the impacts of these factors more clearly, we 
calculate the f 2 of pi in three different ways before mak- 
ing the ratio with the W2 of Pq° ' " : 1) original multi- 



plicity and size, p^°^] 2) multiplicity is adjusted to match 
dN/dy or Pq ° ' ^ , p^° ' " ; 3) both multiplicity and size 

J. i J i .1 Rot,Mul Rot, Mul, Resize mi li 

adjusted to match Pg i Pi . ihe results 

are shown in Fig. [TT]for m — 1 (top- left panel), m — 2 
(top-right panel), m = 3 (bottom-left panel), and m = 4 
(bottom- right panel). Again the k is readjusted inde- 
pendently for each case. We see that matching the mul- 
tiplicity dependence mostly increases V2 at -/Vpart < 200 
where the dN/dy per participant is changing the fastest, 
but has very little influence at TVpait > 200. However, 
after the RMS size of the matter profile is readjusted to 
match that of the participant profile, the calculated V2 
ratios follow the ratio of the eccentricities nicely (except 
in central and peripheral bins for m > 1). We clearly see 
a large sensitivity of V2 on (7^: The large (~ 50%) sup- 
pression of V2 in central collisions, that is, the difference 
between the dashed line and the solid circles, is mostly 
attributable to a '-^ 15% narrowing of ar (Fig. ^. 



o 1 



CO 




0.5 



^^(pRol,»lul,R..iz.y„ (pFlol,»lul 

Ratio of 





100 200 300 100 200 300 







- 


0.2 


r Rot 


'■ 


0.15 


- 


■'^-->% 


J« - 


CM 

> 

0.1 


-ii/ 


/^^^ 




- 


0.05 


p^^ 


^^ 




N, 



part 



N, 



part 



100 200 300 



100 200 300 



N. 



part 



FIG. 11: (Color online) Ratios of calculated V2 between 
Pq°'' " and three cases of Pi^°*: original multiplicity and 
size (dotted lines), multiplicity is adjusted to match dN/dy 
orpo°'' " (dashed lines), both multiplicity and size adjusted 
to match Pq°*' " (solid lines). They are presented separately 
for I\ (top left panel), /2(top right panel), /s (bottom left 
panel) and /4 (bottom right panel). 



cross each other at around A^part ^100. In most central 
collisions the V2 from two component model is suppressed 
by about 10%, even though the eccentricity value shows 
~ 10% increase. Consequently, it under-predicts the data 
in central collisions (see left panel). 



FIG. 10: (Color online) Left: V2 calculated for two Glauber 
geometries, pj^°* and Pi^°*, in their respective rotated frames 
for I\ — 1 4, (from bottom up). Right: Corresponding ratios 
for U2S from I\ — I4, (lines) and for the eccentricity (open 
circles). Note that p^°* (collision density profile) has larger 
eccentricity. 

It is now straightforward to apply what we learned 
from Fig, [m to study the V2 for the two component mat- 
ter profile, p2- P2 is built as a linear combination of par- 
ticipant density profile (pi) and collision density profile 
(po). It is a quite popular initial geometry used in many 
hydrodynamic model calculations [30|, |39|, |40| . It has the 
correct multiplicity; but a smaller geometrical size rela- 
tive to Po due to a centrality-dependent contribution from 
binary collision density profile (Even though 5 is fixed at 
0.14, pcoii becomes more important in central collisions 
because iVcon grows faster than A^part toward central col- FIG. 12: (Color online) Left: V2 calculated for two Glauber 
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lisions.). We naturally expect the corresponding eccen- 
tricity and V2 should sit between that for po and pi . This 
is indeed the case as shown by Fig. [TH The ratio of the 
calculated V2 between p2 and po has shape similar to that 
of the middle curve in Fig. [TI] albeit the rate of change 
is reduced. The V2 ratios decrease with A'part, while the 
ratio of eccentricity increases slightly with A"part- They 



geometries, pj^°' 



and P2^°', in their respective rotated frames 
for I\ —I4 (from bottom up). Right: Corresponding ratios for 
W2S from Ji — 14 (lines) and for the eccentricity (open circles). 
Note that p2 °' (two component model density profile) has 
larger eccentricity. 



As a final note, we point out that the mixing fraction. 



6 = 0.14, between A^part and -^cou is chosen to match the 
dN/dy. However, it is not clear that this parametrization 
necessarily reflects the true shape and size of the matter 
profile (see for example [5(|). This concern is especially 
true in viewing its poor agreement with the central data 
for /i — I4 even after including the eccentricity fluctua- 
tion. Such poor agreement is clearly due to the narrow 
profile of the collision component in the two component 
profile. Nevertheless, it seems that by combining the 
hydrodynamic description of the low px V2 and the jet 
quenching description of the high px V2 , one can gain in- 
sight not only on the eccentricity or shape, but also the 
size of the matter profile. 



D. Fluctuations beyond rotation to the participant 

plane 

The study of fluctuation so far only includes the fluc- 
tuation of PP angle of the initial geometry. In principle, 
we should also consider local density fluctuations which 
affect variance along the long and short axes of the ellip- 
soid without changing the orientation of the PP. In other 
words, both the size Cr and the participant eccentricity 
fpart, which are invariant under rotation, can still fluc- 
tuate event to event for flxed impact parameter. These 
fluctuations are large compared to their mean values, as 
shown by Fig. [131 

Estimation of the influence of these additional fluctu- 
ations requires event- by- event calculation of the jet ab- 
sorption, where the nucleons cannot be treated as point- 
like objects. We assume the nucleon has a Gaussian pro- 
file in the transverse plane with a width of tq in the x and 
y directions, corresponding to a nucleon-nucleon overlap 
function, 



t{x,y) 



1 



2TTrl 



(7) 



and a binary collision profile, 

dx'dy' Ta{x - -,y) X (8) 

TA{x+-+x',y + y')t{x',y') 

The event-by-event participant profile is obtained by 
summing over the nucleon profile for all participants. 

For each event, we generate four dijet pairs by sampling 
its collision profile [Eq. |5] , then calculate their absorption 
in corresponding participant density profile. A total of 
3 X 10^ Glauber events are used. The n is chosen such 
that the overall survival rate averaged over all events is 
0.18 for the 0%-5% centrality bin. To check the stability 
of our result against the finite size assumption of the 
nucleons, we varied the tq from 0.2 — 0.4fm, and we also 
assume the nucleon to be disk of constant density with 
a radius of \J a™^^ ji: jl = 0.58fm. It turns out the final 



results are not sensitive to details of the nucleon overlap 
function, except for very peripheral collisions (A^part < 
20) when nucleon size become comparable to the size of 
the ellipsoid. The deviation is even smaller, when average 
collision geometry is used. More detailed discussion on 
this can be found in Appendix [Ul 

Figure [T3] shows the influences of these additional fluc- 
tuations on jet quenching v^- Again, the jet absorption 
strength is tuned independently to match the central 
i?AA data. The main effect of these additional fluctu- 
ations is a small increase of the V2 in central collisions 
and a small decrease in peripheral collisions. The change 
is less than 10% for TVpart > 100. 

We can understand this 10% centrality dependence 
change of the jet quenching v^ based on what was learned 
from previous discussion as follows. First, we notice that 
the eccentricity averaged over many events is not the 
same as the eccentricity of the matter proflle averaged 



over many events, that is, 



^ 



How- 



ever. Fig [27] in Appendix [C] shows that the difference is 
only about 2% and independent of A'^part and cannot ex- 
plain the difference of the vi in Fig. [TH Thus, it must 
be related to the event-by-event fluctuation of the Ot- 
What really matters is the relative size between the jet 
proflle and the matter proflle, R = ct^^'/ctP''''*. Smaller 
i? implies that more jets are produced in the interior of 
the matter profile, thus suffer more energy loss and have 
large vi ; larger i? implies more jets are produced in the 
corona region of the matter profile, thus suffer less energy 
loss and have smaller V2. 

Figure [151 shows the distribution of R for several cen- 
tralities, as well as its mean value (i?) and standard de- 
viation aji as a function of centrality. (i?) is almost con- 
stant as a function of A^part around ^ 0.85. However 
the width of the distribution ct/e is a strong function of 
centrality; it increases from about 2% in most central 
collisions to more than 5% around A^part ~ 50, and the 
distribution becomes asymmetric toward peripheral bins. 
It is true that the initial jet production rate does not 
depend on the fluctuation of R. However, the survival 
probability does; that is, more jets escape the medium 
when R fluctuates to large values while fewer jets escape 
when R fluctuates to small values. Thus this R depen- 
dent survival probability amplifles the upward fluctua- 
tion of i?, leading to a smaller V2. The suppression of the 
V2 is greater in peripheral collisions due to a broader R 
distribution. This explains the falling of the ratio toward 
peripheral collisions in Fig. [TJ] 

Before closing this section, we stress that the effect of 
fluctuation on V2 can be largely attributed to the fluctua- 
tion of the PP angle. The residual effects, arising mainly 
from event-by-event fluctuation of the relative size be- 
tween matter proflle and jet production proflle, are less 
than 10% for A^part > 100. Thus, it seems reasonable to 
use an averaged matter profile and binary collision pro- 
file for jet quenching calculation, supplemented with a 
small centrality dependent correction. This is of prac- 
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FIG. 13: (Color online) Distributions of participant plane angle (left panel), participant eccentricity (middle panel) and RMS 
size (right panel) for po in three centrality selections (0%-5%, 20%-25% and 50%-55%). 
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FIG. 14: (Color online) Left: V2 for participant density pro- 
file (but in rotated frame and scaled to dN/dy), calculated 
either on averaged profile (Po°*' " ) or event-by-event with 
additional fluctuations (Po°*' "' "'^)- Right: Corresponding 
ratios for V2S from I1 — I4 (lines) and for the eccentricity (open 
circles). Note that the two cases have the same eccentricities. 



tical importance, because event-by-event jet quenching 
calculation is either not possible or computationally pro- 
hibitive for many current pQCD models. However, the 
lumpiness for event- by- event geometry implies large fluc- 
tuation of scatter centers along the jet trajectory, which 
may influence the LPM 



51| coherent effect; there is no 



such problem when the event average density proflle is 
used (see Fig. [511 compare the first three panels with 
the bottom-right panel). However, investigation of such 
effects is beyond the scope of this study. 



namics. Early estimation based on ideal hydrodynamics 
and Glauber geometry [S^l shows that the RHIC V2 data 
require tq < 0.6 fm/c when assuming free-streamin g of 
partons at r < tq. However, Luzum and Romatchke |54| 
argue that the large initial eccentricity of CGC geometry 
allows a bigger tq (up to 1.5 fm/c) for free-streaming, 
without destroying the agreement of their calculation 
with experimental data. Note that free-streaming is an 
extreme assumption, since partons always interact with 
each other and build up flow even if the matter is not 
in local thermal equilibrium [55]. This is especially true 
for jet energy loss which does not explicitly require local 
equilibrium. Nevertheless, it is an interesting question to 
ask whether the high pT anisotropy due to jet quenching 
can provide any constraints on tq. 

Current implementations of the pre-equilibrium energy 
loss are different among various pQCD models. The value 
of To typically varies in 0-0.6 fm/c. Some calculations 
assun ie q = 0, while others assume it is constant at r < 
To [J, [Ij, lla, [53 ■ Although both can describe the single 
inclusive suppression, the extracted q at tq can differ by 
as much as a factor of two [58]. In this work, we tried 
three different formalisms to model the pre-equilibrium 
energy loss. 

• Jets propagate freely to To = Iq/c, then radiative 
energy loss and LPM interference are switched on: 



dl {I - lo) 



I 



dl- 

t + '0 



p(i^ + (; + ;o)v) 



(9) 



E. Dependence on the thermalization time 

One of the main uncertainties in hydrodynamic de- 
scription of the elliptic flow arises from modeling of the 
thermalization time To, i.e. the time at which the system 
reaches local equilibrium and hydrodynamic expansion 
is turned on. This time also explicitly enters the energy 
loss calculation. The value of tq is not known due to lack 
of constraints on the initial geometry and early time dy- 



• LPM effects start at t = 0, but its contribution to 
energy loss is truncated at t < To: 



h = 



I t > to 



/•OO 

/ dl p{-t + B) 

Jla 



(10) 



This functional form is motivated by often made 
claims similar to one from Ref. Il3l: "For times 
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FIG. 15: (Color online) Left: Distributions of event-by-event ratio of RMS width of collision density profile to participant 
density profile, R = CT,':°"/aJ?^'*. Middle: centrality dependence of the mean value for R, (R). Right: centrality dependence of 
the RMS width for R, or. 



prior to 0.6 fm/c, i.e., the .starting point of the RFD 
simulation, we neglect any medium effects, i.e., as- 
sume q = 0. Note that for a purely radiative en- 
ergy loss model where the average energy loss grows 
quadratically with path length in a constant medium 
the effect of initial time dynamics is systematically 
suppressed and no strong dependence of the energy 
loss on variations of the initial time is observed. " 

Radiative energy loss is on all the time, but the 
experienced density is assumed to increase linearly 
with time and reach the local density at tq. 



Ir = 



^dl x/^(^ + '"^ 



\ p{~f + Iv 



) l>lo 



(11) 



This effectively implies a q that grows and reach 
maximum at tq. 

Figure [TC] shows the calculation for various values of 
To from to 2.0 fm/c for the three cases in the rotated 
frame. We find that la exhibits the strongest dependence 
on Tq. The calculated V2 increases almost linearly with 
To and reaches the experimental data at tq ^1.5 fm/c. 
This time is somewhat smaller than a similar analysis 
from Pantuev ^595 , who need To ~ 2 — 3 fm/c to match 
the data, tq is smaller in our case because we include 
eccentricity fluctuation. The increase of V2 with To can 
be attributed to increasingly larger contribution from jets 
originated from the corona region of the overlap, whose 
size is proportional to To [53] ^. Note that the dependence 



^ This can be qualitatively understood as the following. Equa- 
tion |9] defines a corona region I < Iq {la is less suppressed) and 
a core region I > ^o (^a is more suppressed), where I is the dis- 
tance from the surface. The radial distribution of the initial po- 
sitions for the surviving jets is largely defined by the requirement 
that Raa = 0.18 or 18% jets survive in most central collisions. 
When ^0 is small, the corona volume is < 18%, the surviving jets 
come from both the corona and the core and have a quite broad 
radial distribution. As Iq or corona volume grows, more and 
more surviving jets originate from the corona region with larger 



of jet quenching V2 on To is just the opposite of that for 
elliptic flow; the latter always decreases with increasing 
Tq. 

For the second functional form, lb, wc find that the 
truncation of contribution at t < To does simulate the 
suppression of the early contribution due to quadratic 
path length dependence. However, It exhibits a much 
weaker dependence on tq than la, so the two are not 
equivalent. The increase of V2 reaches about 20% for tq = 
0.6 fm/c and grows continuously thereafter. It reaches 
the experimental value at To = 2.0 fm/c instead of To — 
1.5 fm/c for /q. 

As a more realistic scenario that takes into account 
some contributions from t < tq, Ic exhibits much weaker 
dependence on To. The change in V2 is less than 10% at 
To < 0.6 fm/c; but it again increases quickly at large To 
(> 1 fm/c), where matter integral is dominated by the 
corona region which has a large asymmetry. 

A slightly different exploration of effects of early time 
energy loss has been reported in Ref. [58j] . It assumes ei- 
ther a constant q [q = ^(to)] or a g that decreases rapidly 

to its value at to [ (J = ^(to) (^) ]. Both cases imply 
more pre-equilibrium contribution than Ic] thus, we ex- 
pect that they have even weaker dependence on To than 
on Ic. Note that it is generally true that calculation 
which has a smaller to or takes into account the contri- 
bution at T < To always has smaller V2 , because the early 
part of the matter integral tends to be more isotropic. 

Our calculation does not take into account the trans- 
verse expansion. As pointed out earlier [60, [SlI , the de- 
pendence on To is further suppressed if the radial fiow is 
included. This is because the medium moves outward at 
speed of vt- Jets that are generated behind the fluid cell 
need to move and catch up with it. Thus, matter inte- 
grals decrease more slowly with To than the l-flD case. 
Effectively, the radial flow tends to shrink the black core 



anisotropy, and the radial distribution narrows. Until corona vol- 
ume reaches about 18%, most surviving jets come from corona 
and the core becomes almost black. 
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FIG. 16: (Color online) V2 calculated using three different 
schemes for jet absorptions; la according to Eq.|5](top panel), 
lb according to Eq. [10] (middle panel) , I^ according to Eq. [11] 
(bottom panel). In each cases, 8 different thermalization 
times range from /o=0 to 2 fm are shown (curves that are 
lower and darker are for smaller lo). 



region and reduce the dependence on tq. 



F. How sensitive is Raa to initial geometry and 
energy loss formula? 

The setup of our model framework also provides a con- 
venient way to study the centrality dependence of several 
other jet quenching observables, such as inclusive single 
hadron suppression Raa, inclusive away-side per-trigger 
yield suppression Iaa and its associated anisotropy W2** 
(i.e. Iaa as function of angle with respect to the PP). It 
is rather straightforward for us to identify (similar to V2) 
the most relevant control factors of the collision geom- 



etry and path length dependence for these observables. 
To preserve the flow of the main discussion, we focus this 
section on the inclusive Raa, because it coupled directly 
to the V2 discussion. We refer the reader to some initial 
work on Iaa and ^2** i^ Appendix [X] 
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FIG. 17: (Color online) Each of the six panels show the com- 
parison of centrality dependence of Raa between two matter 
profiles (matter types are indicated) for Ji — I4, (eight curves 
in total). The corresponding «2 figure for the same set of 
geometries is indicated in each panel. The first four panels 
shows the effects of switching on and off particular effects of 
the geometry, i.e. Rotation to participant plane for partici- 
pant profile (panel a)), matching dN/dy for participant profile 
(panel b)), matching the size for collisional profile (panel c)), 
and including the additional fluctuation (panel d)). The re- 
maining two panels show comparison between Glauber and 
CGC geometry (Panel e)) and between two Glauber geome- 
tries (Panel f)), respectively. 

Figure [T7] shows the calculated Raa in six different 
cases, with each one designed to check the sensitivity on 
one aspect of the geometry. From top to bottom and left 
to right, the lists of checked effects as follows: 

a) Without {po°^) and with (Po°*' ") eccentricity 
fluctuation. The corresponding V2 comparisons are 
shown in Fig. 3. 

b) Without (p^°') and with (p|^°*''^"') matching the 
multiplicity. The corresponding V2 comparisons are 
shown in Fig. 4. 

c) Without (pf°t^M"i) and with (pR°t^Mui3csize^ ^^^^_ 
justing the RMS size. The corresponding V2 are 
indicated by the dashed and solid lines in Fig. 11. 
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d) Without {pi°*'' ") and with (p^ °*' "' "'^) addi- 
tional fluctuation beyond rotation of the PP. The 
corresponding V2 comparisons are shown in Fig. 15. 



e) Glauber geometry based on (pg °'' " ) versus CGC 
geometry (p3 ). The corresponding V2 compar- 
isons are shown in Fig. 7. 



f) Two glauber geometries: participant profile 
{Pq° ' " ) versus two component profile (/3^°*). 
The corresponding V2 comparisons are shown in 
Fig. 12. 

In each case, all four path-length dependencies /i — £4 
are shown with and without the particular effect under 
investigation. Figure [TH] compares the different choices 
of the thermalization time tq in a broad range using for- 
mulation la (left) and Ic (right); the corresponding V2 
comparisons are shown in Fig. 1161 
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FIG. 18: (Color online) The centrality dependence of -Raa as 
a function of thermalization time lo/c. Two different jet ab- 
sorption schemes are used, that is, la in the left panel (Eq. O 
corresponding V2 figure is the top panel of Fig. I16|) and Ic 
in the right panel (Eq. Illl corresponding V2 figure is the top 
panel of Fig. I16p . Curves that are lower and darker corre- 
sponds to smaller lo- 



It appears that when jet absorption strength is tuned 
to reproduce the suppression in most central collisions, 
the centrality dependence of Raa has limited sensitiv- 
ity on different choices of the collision geometry and 
energy-loss formula. This confirms previous observa- 
tions [13, [531 that the centrality dependence of the Raa 
has limited discriminating power to dynamics of the un- 
derlying energy-loss mechanisms. This is partly due to 
the energy loss bias, but it is also related to the fact that 
-Raa has to vary monotonically between Raa ^ 0.2 at 
large iVpart and Raa ^ 1 when A^part — >■ 0. However, 
Raa does exhibit some sensitivities for three cases; that 
is, it increases somewhat for larger tq (Fig. llSp . stronger 
centrality dependence of dN/dy (Fig. fTTb ) and larger 



RMS size of the matter profile (Fig. [TTb ) ^. However, 
the change is well within typical range of the experimen- 
tal systematic errors. The situation is very different for 
V2 or anisotropy of Raa, which is not a monotonic func- 
tion of iVpart , and exhibits a much greater sensitivity to 
the variation of geometry and energy-loss scheme. 



IV. DISCUSSION AND SUMMARY 

Table U summarizes the sensitivity of f 2 and Raa on 
the choices of matter profile and energy- loss formula (also 
on experimental uncertainties; see Appendix |B]). The 
former includes fiuctuation and shape distortion, which 
affect the eccentricity, centrality dependence of multiplic- 
ity, and the size of the matter profile. The latter includes 
different choices of path length dependence and thermal- 
ization time. Overall, Raa is not sensitive to these fac- 
tors, while V2 has fairly strong dependence on both initial 
geometry and energy loss formula. We can summarize 
the main findings for V2 as follows. 

• Eccentricity fiuctuation has to be included in the 
jet quenching calculation. Without it one can not 
account for the large V2 observed in central colli- 
sions. We estimate that it leads to 15%-20% in- 
crease of the V2 in the mid-centrality bin (Fig. [3]). 

• The residual event-by-event fluctuation, other than 
from eccentricity fluctuation, at most leads to 4%- 
8% change in V2 at iVpart > 100 (Fig.[Tll). 

• A reasonable variation of the multiplicity depen- 
dence, only significantly changes the V2 at A^part < 
100 (Fig.Hl). 

• CGC geometry always results in a larger V2 than 
Glauber geometry at A^part > 150. Depending 
on the choices of Glauber geometry, the increase 
ranges anywhere from 10% (ppart, Fig.[71), 30% (two 
component profile p2, Fig. [T2]), or 20% to > 50% 
(Pcoii, Fig.in]). 

• V2 is very sensitive to change of the RMS width 
(Tr- of the profile. A 15% change in ar can lead to 
about 30%-40% change in calculated V2 fFig. fTTj) . 

These same conclusions seem also to apply for the 
anisotropy of away-side suppression Wj** (see Fig. I^Tj) . 
except that the sensitives seem to be much stronger. We 



^ Naively, one would expect the -Raa to decrease for larger ar be- 
cause more jets originate from inside the profile. However, this is 
only true if jet absorption strength k remains the same. Because 
we always readjust k such that -Raa ~ 0.18 in the most central 
bin, it almost cancels the expected suppression in large -/Vpart 
and even makes the peripheral bin less suppressed (Fig. I17t ). 
This is very different from V2, which always increases for larger 
(Jr (see Fig. llll l. 
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TABLE I: The sensitivity of the V2 and -Raa for various changes in matter profile (relative to participant profile) and energy 
loss schemes. 



Types of changes 


V2 


-Raa 


Fluctuation of RP angle 


Increase by 15%-30% for mid-centrality, much larger 
in central and peripheral (Fig. [3]) 


<5% 


Additional fluctuation 
(relevant for e-b-e jet 
quenching calculation) 


±5% -10% (Fig.lMD 


<5% 


Change in average shape 


Increase by 15%-30% for CGC (Fig. [7] and [HI) , by 10% 
for two component model (Fig. [T2l), 10%-40% for pcoii 
profile (Fig. E]) 


<5% 


RMS size of matter profile 


--0% to -10% change for CGC (compare Fig. [7] to HI), 
^10% to -40% change for pcoii (Fig. |lll compare 
dashed to sohd line) 


10% change for 10% 
difference 


Centrality dependence of 
multiplicity 


Important for A^part < 150 (Fig.g]) 


10%-20% for 
mid-centrality 


Power n for "/" dependence 


Very sensitive 


<5% 


Thermalization time tq 


Sensitive but depends on modeling of energy loss at 
T < tq. For To = 0.6fm/c, it increases by ~ 30% for 
la, < 15% for lb and < 5% for I^. (see Fig.fTB]) 


No, except for 
Apart < 100 


12% variation for Raa in 
0%-5% central bin 


< 7% at iVpart > 100 (Fig. |24|) 


<12% 



stress that the dependence of jet quenching V2 on geom- 
etry is different from that for the low pT V2 driven by 
collective expansion. The latter is sensitive only to the 
eccentricity of matter profile (items 1 and 3 in Table |I| . 

Despite the rather complicated dependence on the ini- 
tial geometry for jet quenching U2, most of them, such as 
the fluctuation and distortion due to saturation, can be 
constrained independently by elliptic flow data [30, [5J] . 
Jet quenching V2 appears to be rather sensitive to the 
choices of energy-loss formula, thus making it an ideal ob- 
servable for gaining insights on energy-loss mechanisms. 
Based on the comparisons shown in Fig.[3|and[7J it seems 
that naive path length dependence motivated by radia- 
tive energy loss, /i, is insufficient to describe the data for 
both Glauber and CGC geometry even with eccentricity 
fluctuations taken into account. It appears that either 
2 < TO < 3 for Glauber geometry based on participant 
profile or TO ~ 2 for CGC geometry, both with eccen- 
tricity fluctuations, have the best match with the data. 
Note that m = 2 corresponds to the AdS/CFT type of 
energy loss AE ex l^ for a strongly coupled plasma [23 |. 
Similar strong path length dependence is also observed 
for away-side suppression Iaa and it's anisotropous V2** 
(seeFig.[iniand[2I|). 

The jet quenching V2 is also quite sensitive to thermal- 
ization time To, but the sensitivity depends on the model- 



ing of the pre-equilibrium energy loss. By assuming free- 
streaming up to To — 1.5 fm/c and including eccentricity 
fluctuation, we can reproduce the experimental data with 
quadratic path- length dependence of energy loss. How- 
ever, inclusion of a very modest pre-equilibrium energy 
loss, for example, Ic which assumes a q which linearly 
grow to ^(to) at To, already signiflcantly suppressed the 
dependence on tq up to 1 fm/c. We can draw similar con- 
clusions for away-side suppression Iaa and the associated 

,,-fAA 



anisotropy W2** (see Figs. [Hj and 

The preceding discussion attested to the value of jet 
quenching V2 in understanding the roles of various geom- 
etry factors and constraining the energy loss mechanisms. 
One can obtain more discriminating power by combining 
all jet quenching observables, V2, Raa, Iaa, and V2^^ 
(see Appendix [Xj for more discussions) . Initial theoret- 
ical work already demonstrated the value of combining 
the i?AA and Iaa [8*,^]; One can do a better job by also 
including a calculation of V2- It is worth pointing out 
that the study of high px V2 benefited significantly from 
extensive experimental and theoretical work on the low 
Pt V2, which has provided important constraints on the 
initial eccentricity. We have shown that the reverse is 
also true; i.e. jet quenching observables can also provide 
useful new insights on the initial geometry, which could 
help the interpretation of the elliptic flow data. 
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In summary, using a simple jet absorption framework, 
we studied the sensitivity of jet quenching V2 to various 
aspects of colhsion geometry and the path-length depen- 
dence of energy loss. Besides the eccentricity, we found 
two other ingredients of the collision geometry, namely, 
the centrality dependence of the matter integral and the 
relative size between the matter and the jet profiles, are 
important for jet quenching V2- We compare the cal- 
culated V2 from both Glauber and CGC geometry with 
experimental data. A path length dependence stronger 
than the native AE ex l^ dependence from radiative en- 
ergy loss, as well as the inclusion of the eccentricity fluc- 
tuation, are necessary to reproduce the V2 data. A de- 
tailed comparison between Glauber and CGC geometry 
shows that a 15%-30% increase of initial eccentricity in 
CGC only results in half the increase in calculated V2 
due to a small narrowing of the CGC geometry. This 
points to an interesting possibility: A large V2 can be 
easily generated if the jet production proflle is signifi- 
cantly narrower than the matter profile. This requires 
a transverse profile distribution that narrows with mo- 
mentum or Bjorken variable x = 2pt/\/s- This happens 



for the nucleon-nucleon collisions |62l |63[; however, we 
are not aware yet of a physical mechanism to produce a 
significant narrowing at large x in heavy nuclei. 

Our estimations are based on a simple jet absorption 
framework. Admittedly, it is too simplistic to give di- 
rect insight on the dynamics of the energy-loss process. 
However, it proves to be a useful tool for understanding 
the centrality dependence of various jet quenching ob- 
servables, for identifying the most relevant factors in the 
collision geometry and path-length dependence, as well 
as for estimating the sign and magnitude of the change 
as we vary those factors. We have documented all the 
matter and jet profiles used in this study in Ref. [6J]. 
They can be used as input for future, more-realistic jet 
quenching calculations. 
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Appendix A: Jaa and its azimuthal anisotropy? 

Initial geometry should also leave footprints on the az- 
imuthal distribution of the away-side jets. The appropri- 
ate observable for this purpose is the anisotropy of the 
per-trigger yield for the away-side jets, /AA(0trig-*part), 
which refiects the path length dependence of the energy 
loss for the away-side jet. Due to the surface bias of the 



trigger jets, the away-side jets on average have longer 
path length to traverse; thus, they are expected to ex- 
hibit stronger suppression and larger anisotropy. In this 
work, we calculate the anisotropy coefficient as 



/a A 



(/AACOs2(0trig-*part)) 



(Al) 



As a warm-up exercise, Fig. llQI shows the azimuthal de- 
pendence of per-trigger yield suppression /aa in 0%-20% 
and 20%-60% centrality, calculated for Po°*' " profile, 
i.e. participant density profile in rotated frame and re- 
scaled to match the dN/dy data. As one increase the 
power m of the path length dependence in Im, Iaa for 
the 20%-60% bin shows a dramatic decrease in the out-of- 
plane direction where the path length is large, but only 
a modest decrease in the in-plane direction where the 
path length is small. That is because the suppression 
for a given centrality is largely determined by the typical 
matter integral Im ~ (L)"^ / (i™)™, which changes more 
rapidly for larger m. Here L™ is some typical length scale 
fixed in central collisions for m. Clearly, the increased 
sensitivity of large m can generate a large anisotropy, 
hence large U2**- However, the price one has to pay is 
that it leads to a large suppression in the 0%-20% bin 
as shown in the left panel of Fig. [191 This is because a 
large W2** naturally implies a strong suppression in cen- 
tral collisions, as long as the suppression is a monotonic 
function of the path length. 
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FIG. 19: (Color online) The lAA{(ptrig - *part) for 0%- 
20% (left panel) and 20%-60% (right panel) centrality bins. 
The points are results of calculation, the line is a fit to 
1 -I- 2ii2** cos 2A4> + 2^4** cos 2A<j>. We see a finite V4** sig- 
nal, but it is less than 10% of 112**, and is ignored in this 
study. 



The jet absorption framework used in this work is a 
pure geometrical model in that the calculated suppres- 
sion depends only on the path length, thus it always pre- 
dicts Iaa < Raa due to the longer path length of the 
away-side jet. However, experimental data seem to sug- 
gest that Iaa ^ Raa 0, Isl, [6a] . This is because that 
the away-side associated hadron spectra are much flatter 
than the inclusive distribution, due to the requirement 
of a high pt trigger. For a typical trigger of 5 GeV/c, 
the away-side conditional spectra in p + p collisions at 
yG = 200 GeV, if parameterized via a power law func- 
tion 1/pJ, have a power of n = 8, in contrast to n = 4 for 
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inclusive hadrons [63|- A simple estimation shows that 
to reach the same level of suppression, it takes about 
50% more energy loss for away-side jets than for inclu- 
sive jets '67]. What this means is that treating energy 
loss as absorption is not sufficient; an energy shift term 
is required as well. However, phenomenologically, we can 
still use the jet absorption framework, if we allow the k 
to also depend on power n; i.e, a smaller k is required for 
an away-side jet due to a flatter input spectra. Such de- 
pendence in principle can be fixed by the /aa data from 
STAR and PHENIX g,!!!^. However, we shah defer 
this improvement to a future study. 

Nevertheless, the current setup is sufficient for study- 
ing the sensitivity of Jaa(A^) on the choices of under- 
lying collision geometry. Figures [501 and [5T] summarize 
the relative change of the inclusive /aa and its anisotropy 
V2^^ , respectively, as one varies various aspects of the col- 
lision geometry. The six cases in both figures correspond 
to the same change in collision geometry as in Fig. [T7] 
for inclusive Raa, that is, a) eccentricity fluctuation, b) 
matching the multiplicity, c) re-scaling the geometrical 
size CTr, d) additional fiuctuation not included by rota- 
tion, f ) default Glauber geometry versus CGC geometry, 
g) default Glauber geometry versus two-component ge- 
ometry. 

In general, we see that /aa and ^2** are much more 
sensitive to their counterpart single particle observables, 
Raa and V2- They are also very sensitive to both the 
path-length dependence and the thermalization time tq 

,-fAA 



(see Figs. [22] and 



In one case (Fig. [2Tb)). the 



even becomes negative in central collisions, reflecting the 
dominance of tangential emission when the collision pro- 
file is used as the matter profile (pi). This is the case 
because pi has the narrowest profile (15% smaller than 
participant profile), such that more jets are generated 
on the surface and can survive if emitted tangentially. 
Unlike the V2, the change in ^2**, when tuned to same 
multiplicity and average, is not proportional to the corre- 
sponding change in e. The only change to which /aa and 



.Jaa 



is not sensitive is the fluctuation of the PP angle. 



Appendix B: Sensitivity to k 
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FIG. 20: (Color online) Conventions are similar to those in 
Fig. 1171 Each of the six panels shows the comparison of cen- 
trality dependence of /aa between two matter profiles (mat- 
ter types are indicated) for /i — I4, (eight curves in total); 
the corresponding V2 figures for the same set of geometries 
are indicated in each panel. The first four panels shows the 
effects of switching on and oif particular effects of the geom- 
etry, i.e. rotation of participant plane for participant profile 
(panel a)), matching dN/dy for participant profile (panel b)), 
matching the size for coUisional profile (panel c)), and includ- 
ing the additional fluctuation (panel d)). The remaining two 
panels show comparison for Default Glauber vs. CGC geom- 
etry (Panel e)) and between two Glauber geometries (Panel 
f)), respectively. Note that the small difference in panel a) 
between with and without rotation for m = 4 should be at- 
tributed to increased sensitivity of /aa to changes in -Raa for 
large m (because n is tuned by hand to match Raa ~ 0.18 
only to the third digit after 0). 



We have discussed four physics observables, Raa, ^2, 
/aa, and W2**- Their centrality dependencies and path 
length dependencies are entirely controlled by the k pa- 
rameter introduced in Sec. [Hj For each running mode, 
K is always readjusted independently such that Raa = 
0.182 in the 0%-5% centrality bin. 

The discussion so far has ignored the experimental un- 
certainty on the Raa measurement, and it is about 12% 
for the 0%-5% centrality bin (0.160 < Raa < 0.204). 
This uncertainty invariability translates into some uncer- 
tainty on K, which in turn translates into some uncertain- 
ties on the other three observables. Table. [Ill summarizes 
the K values and associated uncertainties calculated for 
po, p2 and ps and m = 1 — 4. One can see that the 



fractional uncertainties increase gradually for larger m. 
The differences between po, P2 and ps, however, can be 
largely attributed to their different total matter integral 
in the 0%-5% bin, which are 353, 301 and 695 for po, p^ 
and Pa, respectively. If we readjust the total matter in- 
tegral for the 0%-5% centrality bin to the same number, 
say 353, then the obtained k would need to be multiplied 
by 301/353 = 0.854 for p2 and 695/353 = 1.97 for pa, 
respectively. This procedure makes the n rather close to 
each other for different p's; the residual differences then 
can be attributed to their different shapes and sizes. 

Figure [M] shows the centrality dependence of the four 
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FIG. 22: (Color online) The centrality dependence of /aa as 
a function of thermalization time lo; Two different jet absorp- 
tion schemes are used, i.e. la in the left panel (Eq. [9| and 
Ic in the right panel (Eq. Ill|) . Curves that are higher and 
darker corresponds to smaller lo (opposite to the ordering for 
V2, Raa and ^2**)- 
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FIG. 23: Similar to Fig. 1221 except plotted for v. 
that are lower and darker corresponds to smaller lo 



Curves 



jet quenching observables using a po niatter profile ^ with 
K varied to allow the Rkk to change by ±12% in the 0%- 
5% centrality bin for each vn value (panel (a)). These k 
values are then used to predict the centrality dependence 



FIG. 21: (Color online) Conventions are similar to those in 
Fig. 1171 Each of the six rows show: (left panel) the com- 
parison of centrality dependence of v^^ between two matter 
profiles (matter types are indicated), (right panel) The cor- 
responding ratios together with ratios of eccentricities. The 
corresponding ^2** Plot to'" the same set of geometries is ndi- 
cated in each panel. The first four rows show the effects of 
switching on and off particular effects of the geometry, i.e. Ro- 
tation to participant plane for participant profile (Row a)), 
matching dN/dy for participant profile (Row b)), matching 
the size for coUisional profile (Row c)), and including the ad- 
ditional fluctuation (Row d)). The remaining two rows show 
comparison for Default Glauber vs. CGC geometry (Row e)) 
and between two Glauber geometries (Row f)), respectively. 



TABLE II: The k values for three matter profiles, po (partic- 
ipant profile), p2 (two component profile with 5 — 0.14), and 
P3 (CGC profile), and four / dependencies (m = 1,2,3,4). 
The range of k in each case corresponds to 0.160 < -Raa < 
0.204 in the 0%-5% centrality bin. The total matter integrals 
for each profile are also listed. 



Po 

P2 

p-i 



Integral 



353 
301 
695 



h 



0.147 



-10% 



0.185^ 



+ 14% 
11% 

0.076+^g^° 



0-082+i«;: 
117+20% 

0.046+1^^: 



003511^ 
0.054+?e^: 
0.020+^^^° 



0.0125+^2?% 
0.0206+^22% 
0.0073+2?% 



"■ Similar results are also observed for p2 and pa matter profiles. 
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of V2 (panel (b)), /aa (panel (c)) and fj** (panel (d)). 
We can see that the fractional change on V2 is typically 
smaller than 5% for A^part > 100, and the change de- 
creases for larger A'part; the direction of the change, how- 
ever, is anti-correlated with the change of _Raa- The frac- 
tional change on /aa is about 20%, about twice amount 
of the change of Raa\ this change is in the same direc- 
tion as that for the i?AA- The fractional change in U2*'^ is 
more complicated: It typically decreases for central col- 
lisions and increases for mid-central and peripheral colli- 
sions, but the overall change is typically less than 5%. 
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FIG. 24: (Color online) Centrality dependence of (a): Raa, 
(b) ii2, (c) /aa and (d) ^2** using a po matter profile with 
K varied to allow the Raa to change by ±12% in the 0%- 
5% centrality bin. Different line styles are used for m — 1 
(solid lines), m = 2 (dotted lines), m = 3 (dashed lines), 
and m — 4 (long dashed lines); the upper (lower) values of k 
are represented by thin (thick) lines. Similar results are also 
observed for p2 and pa matter profiles (not shown). 



Appendix C: Comments on finite nucleon size effect, 
event-by-event fluctuation etc 

We can show explicitly why the finite nucleon size is 
not important for V2 calculation except in most periph- 
eral collisions. We notice that the finite nucleon size leads 
to an increase of the variance matter profile, but does not 
change the orientation of the rotated frame: 



The eccentricity decreases a little bit due to smearing of 
n-n overlap function. However, for a typical fireball size 
of ar = 3fm, this is only a 1.7% change in the eccentric- 
ity. ^ The RMS size of the ellipsoid also increases by a 
few percent, but the corresponding jet production pro- 
file increases by the same amount, resulting in almost no 
change on the calculated V2 ■ Figure [25] shows the ratio 
of the V2 for event averaged po calculated assuming the 
nucleon profile follows either the 6 function or a Gauss 
function with width of tq = 0.4 fm. The differences of 
V2 are well within 3%, except at TVpart < 20, where the 
To = 0.4 fm case is smaller. 
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FIG. 25: (Color online) Comparison of the t;2's for Glauber 
geometry, which are filled either according to center positions 
of nucleons or assuming Gauss profiles with a width of 0.4 fm. 

Note that similar findings also appear in MC-KLN, 
which assumes that participating nucleons are disks of 
finite size in filling the participant profile distribution (in 
addition to assuming finite nucleon size for determining 
whether it is a participant). This is illustrated by Fig. [511 
which compares the V2 results for the same Glauber ge- 
ometry determined with either the PHOBOS code or the 
MC-KLN code. The ratio drops at iVpart < 20 similar to 
Fig. [23 As a side note, we point out that this plot also 
confirmed the consistency between the PHOBOS code 
and the MC-KLN code for calculating the Glauber ge- 
ometry when running with the same parameters. 

We have seen that the event-by-event fluctuation leads 
to large dispersion of the distribution of various geometri- 
cal variables, such as eccentricity and RMS size in Fig.lT^ 
Yet the mean values seem to be insensitive to whether 
they are calculated for each event and then averaged over 
many events or are calculated directly from the averaged 
profile. As an example, Fig.[57|compares the eccentricity 

averaged over values calculated event by event ( Htt — ? 

with the eccentricity calculated from the average profile 

jprr^TpT)- The ratios are shown in the right panel for 



= 

^'2 



' xy 
^'2 



The impact is somewhat larger for erp, i.e. average without 
rotation 
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both eRp and Cpart- As one can see, the difference is < 2% 
for Cpart and somewhat larger for e^p ^. This is quite 
remarkable given the sharp visual contrast between the 
lumpiness of the event-by-event profile and the smooth- 
ness of the average profile (see Fig. 
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Nevertheless, the good news is that the average profile 
seems to preserve most of the relevant geometrical infor- 
mation. For example, one can calculate the eccentricity 
directly from the overall matter profile instead of calcu- 
lating it event by event and then averaging it over many 
events. 



FIG. 26: (Color online) Comparison of the t;2's for Glauber 
geometry calculated using the PHOBOS code and the MC- 
KLN code. 
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FIG. 27; (Color online) Left panel: the standard (circles) and 
participant (boxes) eccentricity either by averaging over their 
event-by-event value (filled symbols) or calculated directly 
from the averaged matter profiles (open symbols). Right 
panel: the corresponding ratios for the standard and par- 
ticipant eccentricities. 

"^ The agreement generally worsens when the width of the distri- 
bution becomes large relative to the mean value, which is the 



case for CRp. 
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rotated to participant plane. 
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